Three-Dimensional Phase-Kink State in Thick Stack of Josephson Junctions 

and Terahertz Radiation 
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The dynamics of superconductivity phase in thick stack of Josephson junctions with strong in- 
ductive coupling, such as the one realized in layered high-T c cuprates and possibly the recently 
discovered FeAs-based superconductors, is investigated under a c-axis bias voltage and in the ab- 
sence of an external magnetic field. The kink state found previously by the present authors is 
extended to three dimensions for both rectangular and cylindrical geometries. The IV characteris- 
tics are calculated and the distributions of electromagnetic field inside the samples are clarified. The 
solution for a cylindrical mesa exhibits a higher resonating frequency than that of a square mesa 
with the same linear size by a factor of ~ 2.4. More importantly, from the radius dependence of the 
resonance frequency for the cylinder geometry it is possible to confirm directly the kink state, and 
thus to reveal the mechanism of the strong radiation discovered in recent experiments. 

PACS numbers: 74.50.+r, 74.25.Gz, 85.25.Cp 



I. INTRODUCTION 

Stimulated by the recent progress in radiation of ter- 
ahertz (THz) electromagnetic waves using a mesa of 
Bi 2 Sr 2 CaCu20 8+< 5 (BSCCO) single crystal [10, we in- 
vestigated the dynamics of superconductivity phase of 
the intrinsic Josephson junctions with strong inductive 
couplings. A new dynamics state was found in the ab- 
sence of an external magnetic field, in which the gauge 
invariant phase difference in each junction is character- 
ized by static ±7r phase kinks stacked periodically along 
the c axis, in addition to a phase term evolving linearly 
with time accordingto the ac Josephson relation and 
the plasma term [3j, |j| . The phase kink permits pump- 
ing large dc energy into the plasma oscillation via the ac 
Josephson effect, and enhances strong radiations of THz 
electromagnetic waves. This state is the only one known 
so far which counts for the following important features 
of the coherent radiations observed in Refs. Q, Qj: (1) 
the frequency is determined by the cavity mode of the 
mesa sample, (2) the frequency and voltage satisfy the 
ac Josephson relation, and (3) the radiation takes place 
in a sharp regime of voltage. 

For simplicity a two dimensional (2D) solution was 
worked out explicitly, presuming a state uniform along 
the in-plane direction of the long edge of the mesa j3|. 
This solution corresponds to the (1,0) cavity mode. In 
the present work, we show that the kink state can be 
extended to 3D, where the phase kink runs in the ab 
planes, and the solution is labeled by the (1,1) cavity 
mode. The kink state can also exist in samples of the 
cylinder geometry. 

We first derive a set of equations for the shape of the 
phase kink, the amplitude and phase shift of plasma oscil- 
lation. An approximate solution can be obtained easily in 
the strong inductive- coupling limit where the phase kink 
renders a step function. This formalism releases one from 



a full computer simulation, and can be extended from 2D 
to 3D straightforwardly. 

The present paper is organized as follows. In Sec. II, we 
formulate the kink state for the (1,0) cavity mode (thus 
2D) in a rectangular mesa presented in Ref. [3j . The anal- 
ysis is then extended to 3D in Sec. Ill and the (1,1) cavity 
mode is analyzed. The current-voltage curve of the kink 
state is presented, and the distribution of electromag- 
netic field is derived. Section IV is devoted to analysis 
for the cylinder geometry. Radiation of THz electromag- 
netic waves is addressed under a boundary condition of 
an effective impedance in Sec. V. Discussions are given 
in Sec. VI before a brief summary. 



II. ID KINK IN RECTANGULAR MESA 

The coupled sine-Gordon equations for a stack of 
Josephson junctions p,[a,01 take the following form when 
the inductive coupling C, is large 
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where Pi is the total gauge-invariant phase difference 
of the Z-th junction, the lateral directions are scaled 
by A c , time by the plasma frequency to p , external cur- 
rent by the critical Josephson current J c , C, = A 2 6 /sD, 
j3 = AircT c \ c /cy/r c and A^Qi = Q l+1 + Q ; _i - 2Q h with 
X a b and A c the penetration depths, s (D) the thickness of 
the superconducting (insulating) layer, a c the c-axis con- 
ductance due to quasi particles, e c the dielectric constant, 
and c the light velocity in vacuum [3,|g|. We consider first 
the system under the boundary condition d n Pi = with 
n the normal to the sample edges 0, [i| ■ 

For the kink solution of the (1,0) cavity mode of a 
rectangular mesa, the phase difference Pi is given by [3] 



Pi(x,t) = cut + A cos{—) sm(ut + ip) + fiP s (x), (2) 

where fiP s (x) is a kink from to ±7r with the center at 
x = L x /2. The first linear term is in accordance with the 
ac Josephson relation. 

In a state characterized by the phase Eq. ([2]) with /; = 
(-1)' or fi = (-l)!'/ 2 !, Eqs. H) are decoupled. From 
the static parts, the relation of current conservation is 
derived Q: 
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where g\ = cos(ttx/L x ) and J„ is the Bessel function of 
the first kind. The phase kink P s (x) is governed by the 
differential equation 



d 2 x P s = q( cos (pJ-i(Ag[ ) sinP s 



(4) 



with q = 4 for period-2 and q = 2 for period-4 solu- 
tions, and the boundary condition d x P s = at the edges. 
While this differential equation has to be solved numer- 
ically, to a good approximation, the phase kink can be 
described well by P s (x) = 4>{x)/2 with the soliton solu- 
tion 4>{x) — 4arctan[exp((x — L x /2)/w)}. The width of 
the phase kink is w ~ l/^/Z ~ 10~ 3 for BSCCO, which 
renders the kink almost a step function. 

From terms associated with sin(wi) and cos(ut) when 
the phase Eq. ^ is subject to Eq. ([T]), one obtains 



smyff' j^gl (J (Agl ) + J_ 2 (Ag{ )) cos P s 



\A{^ - (n/L x f) = 
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With the aids of Eq. (O and the recursion relation of the 
Bessel functions zJ v _i(z) + zJ„ +1 (z) = 2vJ v {z), Eq. © 
is simplified as 
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where the first term clearly counts for the dc ohmic dis- 
sipation, and the second term is caused by the ac voltage 
generated by the plasma oscillation, which modifies the 
otherwise linear IV curve. 

The four equations ((4j)-(|Tj) describe the kink state spec- 
ified by five quantities, namely J cx t, w (i.e. voltage), A, 
ip and P s (x). Especially, the IV characteristics can be 
derived by sweeping the voltage u. 

As the frequency, or equivalently the voltage ap- 
proaches the cavity mode u> = n/L x , large plasma os- 
cillations are stimulated and a large external current is 




FIG. 1: (color online). 2D phase kink for the (1,1) cavity 
mode of a rectangular mesa obtained by numerical integration 
on Eq. HID). 



shunted by the Josephson current. Just at the cavity 
frequency Eq. ^ is reduced to 



— / dxgl (J (Agl ) - J- 2 (Ag[ Q )) cos P s = 0. (8) 

J-ix Jo 

In the limit of strong inductive-coupling, where the phase 
kink is well approximated by a step function, Eq. (JSJ) is 
easily solved to yield A — 2.331, and then smip — 0.4625 
and J ext -0co = 0.2133 for (3 = 0.02 and L x = 0.4 in units 
of A c = 200 fxm, typical for the BSCCO single crystal used 
in recent experiments. 

For comparison, simulations gave J cxt — (3u> — 0.37 @]. 
Therefore, the present simplified equations give a semi- 
quantitative evaluation for the kink state at the cavity 
modes. 



III. 2D KINK IN RECTANGULAR MESA 

Now we look for the state of static phase kink of two 
lateral dimensions, which corresponds to the (1,1) cavity 
mode of a rectangular sample: 
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The equation for the 2D kink is 



(dl + 8 2 y )P s = q{ cos ipJ-iiAfu) sin P s , (10) 

where g\ x = cos(irx/L x ) cos(7ry / L y ) with d n P s = at 
edges where n is the normal of the edges. The solution 
obtained by numerical integration is depicted in Fig. [TJ 




FIG. 2: (color online). IV curves derived by Eqs. (fTTj) - ([T3)) 
presuming a step-function phase kink (red symbols) and by 
Eq. (|14p (blue lines) for the (1,1) mode of a rectangular mesa 
with L x = 0.4 and L y — 1.5. The cavity frequency is u — 
8.128, which corresponds to 0.62 THz approximately. 



with L x — 0.4 and L y — 1.5 as in experiments [l|, and 
all the prefactors arc included into the width of the kink. 
The remaining three equations are given by 



\Afiu = 

smtptf* f L " ^Lg^MAgl,) + J- 2 {Mi)) cosP* 



(11) 



4^% 2 -(xr) 2 -(t) 2 ) = 

cos^ / L » gftglMAglJ- J- 2 (49ii))cosP*, 

(12) 




FIG. 3: (color online). Distribution of (a) electric field and (b) 
magnetic field in for the (1,1) cavity mode of the rectangular 
mesa. 
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Approximating cosP s by the 2D step function (see 
Fig. [I]), the IV curve is evaluated and displayed in 
Fig. [21 An enhancement of the input current appears 
when the voltage approaches the cavity mode value w — 
^(tt/L x ) 2 + (ir/Ly) 2 , where the extra energy is pumped 
into plasma oscillation. 

The way of the system evolving into the resonance 
state from the linear ohmic regime can be seen more 
transparently by neglecting A in the right-hand side of 
Eqs. CO} and (T2"]). One then finds 
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with I{ t — (l/L x L y ) L x J Q y dxdygli cosP s . A similar 
expression has been derived for the ID kink state [J] . The 



development of the resonance is captured analytically by 
the above expression, although it becomes invalid around 
the resonating regime where A is not small any more. In 
the strong inductive-coupling limit where the phase kink 
takes a step function, P n = 4/ir 2 . 

As seen in Fig. [21 the two treatments give the same IV 
curves outside the resonance regime where the plasma 
amplitude is small. Usage of Eqs. (fTT]) and (fT2| improves 
the description of the resonating state over Eq. (fT4|) . 

The electromagnetic fields are given by E z (x,y) = 
d t p\x,y), B x (x,y) ~ -d y P(x,y) and E.y(x,y) ~ 

d x P(x,y), with P(x,y) the spatial part of the plasma 
term in Eq. (J5J) 3]. As shown in Fig. the electric 
field takes the maximal absolute value at the four cor- 
ners, whereas the magnetic field is maximal at the cen- 
ters of the edges. The magnetic field penetrates into the 
system from the two edges along the x direction, and 
flows away from the two edges along the y direction, and 
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FIG. 4: (color online). Distribution of static Josephson cur- 
rent for the (1,1) cavity mode of a rectangular mesa. 



FIG. 5: (color online), (a) Phase kink and (b) radial depen- 
dence of phase in a cylinder sample obtained by numerical 
integration on Eq. (|16|l . 



the pattern oscillates with time by switching x and y. 
The y-component of magnetic field is larger in absolute 
value than the ^-component since the system is longer 
in the y direction. Both the electric and magnetic fields 
are shielded in the superconducting sample and thus the 
central part of the mesa is free of electromagnetic fields. 
As discussed in Ref. Q , in the resonance state at cavity 
modes higher harmonics exist and the simple symmetry 
in Fig. is modified. 

The static phase kink generates static Josephson cur- 
rents in junctions as depicted in Fig. |4j The patterns 
with swapped up- and down-ward Josephson currents 
are stacked periodically along the c axis according to 
ft = (-1)' or f, = (-l)P/2]. 



\APuJ 2 (v m ) = 
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(19) 



where Jo(voi) = —0.4028, and the normalization relation 
Jo J(){ v oip/ a )pdp = a 2 J 2 («oi)/2 has been used. The fre- 



IV. KINK STATE IN CYLINDER GEOMETRY 

Now we consider a sample of cylindrical shape and ra- 
dius a. For the (0,1) mode in the cylinder geometry which 
is isotropic azimuthally, the total phase difference is given 
by 



Pl(p,t) =ut + AJ {^) sin(wi + v>) + fiP s (p), (15) 
a 

where vqi — 3.8317 is the first zero of J\. 

The equation for the static phase kink in cylinder is 



(d 2 p + (l/p)d p )P s = q( cos ^J-i(^oi) «n P s , (16) 

with g^ = Jo(voip/a) under the boundary condition 
d p P s = at the edge. The solution obtained by numer- 
ical integration is displayed in Fig. [5] for a = 0.4. The 
phase crosses n/2 at p = u a ia/v i where it i = 2.4048 is 
the first zero of Jo- 

The remaining equations are given by 




FIG. 6: (color online). IV curves derived by Eqs. (fl"7)) - (tT9)) 
presuming a step-function phase kink (red symbols) and by 
Eq. (|20[l (blue lines) for the (0,1) mode of a cylinder sample 
of a = 0.4. The cavity frequency is ui = 9.579, which corre- 
sponds to 0.73 THz approximately. 



quency of the (0,1) cavity mode for the cylinder geometry 
is u) = VQi/a which reaches ITHz when a — 58.5/im. 

The IV curve is displayed in Fig. [6j taking cos P s as 
the step function. For comparison we also show the IV 
curve evaluated by the following expression 



J cxt = /3w(l + 



(w 2 -(« i/a) 2 ) 2 + (/^) 



r). 



(20) 



with /q-l = (1/a 2 ) J Q pdpg^ cosP s obtained in the same 
way as Eq. (fT4|) . For the step- function phase kink, Ifa — 
0.1701. Once again the two evaluations agree with each 
other outside the resonance regime. 

The electromagnetic fields take the following forms: 
E z = AwJ (v ip/a), B x = A{v m y /2pa)Ji{v m p/ a) and 
B y = —A(v n ix/2pa)Ji(voip/a). The distribution of elec- 
tromagnetic field is shown in Fig. [7] which oscillates with 
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FIG. 8: (color online), (a) Distribution and (b) radial depen- 
dence of Josephson current for the (0,1) mode of a cylinder 
sample. 
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time. The electric field takes the maximum at the center 
of cylinder, changes its sign at p = woia/woi = 0.6276a 
and at the perimeter it reaches —0.4028 of the maximum 
value (Figs.EIa) and (b)). The magnetic field is circular, 
the amplitude is zero at both the center and edge of the 
cylinder, and presumes its maximum at ~ a/2 (Figs. [7(c) 
and (d)). The distribution of the static Josephson cur- 
rent generated by the static phase kink is displayed in 
Fig. [5] The patterns with swapped up- and down-ward 
Josephson currents are stacked periodically in the c axis. 



V. RADIATION OF ENERGY 

The boundary condition of a single junction subject to 
radiation was discussed in literatures [9|, [l(J, LUj • There 
is a significant mismatch of the impedance between the 
junctions and outside space. As in the previous work [3j, 
we model the space by an effective impedance Z. The 
effective impedance should be very large due to the small 
thickness of the sample compared with A c [12j . We have 
confirmed that the cavity resonance based on the kink 
state is stable for \Z\ > 50. 

The IV characteristics should be modified when the 
radiation is present. Since the radiation is governed by 
the attached impedance, the IV characteristics can be 
evaluated simply by an effective parallel circuit. With a 
simple dimension counting one finds for a real impedance 

Z 



(3u + /3ujA 2 /A + uA 2 /{ZL x ), 



(21) 



FIG. 7: (color online). Distribution of electric field ((a) and 
(b)) and magnetic field ((c) and (d)) in a cylinder sample. 



J oxt = /3uj + /3uA 2 /8 + loA 2 {1/L x + l/L y )/{2Z), (22) 



VI. DISCUSSIONS AND SUMMARY 



J cxt = 0u> + (3u;A 2 J 2 Ki)/2 + uA 2 J 2 (v i)/(Za), (23) 

for the (1,0) mode uniform in the direction of long edge 
and the (1,1) mode of the rectangular mesa, and for the 
(0,1) mode of the cylindrical mesa, respectively. The 
density of radiation energy measured by the Poynting 
vector then reads 



P 10 = w*A'/(2Z), 



PI, = u 2 A 2 /(4Z), 



PS 1 =co 2 A 2 J 2 (v 01 )/(2Z), 



(24) 



(25) 



(26) 



for the respective modes. 

Presuming the same impedance, we evaluate the ratio 
between the Poynting vectors at the corresponding cavity 
modes when the linear sample size is the same 2a = L(= 
L x = L y ): Pgx/Pfo = 3.125 and Pfo/PJa = 1.898, where 
the maximal amplitudes A = 2.331 for (1,0) mode and 
2.991 for (1,1) mode of the rectangular mesa, and 4.194 
for (0,1) mode of the cylinder geometry have been used. 

There are three factors in determining the Poynting 
vector: the maximal plasma amplitude, the factor for 
the electromagnetic field at the edge, and the cavity 
frequency. In cylinders the electromagnetic fields are 
isotropic azimuthally, and can make a full use of the 
phase kink in pumping dc energy into plasma oscilla- 
tions, which results in a large plasma amplitude. The 
(1,1) mode in the rectangular mesa is also better than 
the (1,0) mode since the former can use a 2D phase kink. 
While the electric field is maximal at the edges of the 
rectangular mesa, the edge electric field for the cylinder 
geometry is only 40% percent of its maximum taken at 
the center of the cylinder. This suppresses the radiation 
power as in Eq. (|26[) . The cavity frequency in the cylin- 
drical geometry is superior to the rectangular geometry 
by a factor ~ 2.4 for the same linear system size. 

While the cylinder geometry exhibits a large density 
of power emission at the sample surface, the energy is 
radiated uniformly in the azimuthal direction. In appli- 
cations one needs to focus to a small angle. In contrast, 
the directivity of the radiation from a rectangular mesa 
is helpful for gathering energy efficiently. 

As can be seen in Eq. ([5]) and the counterparts for 
the other modes, the plasma amplitude at the cavity 
resonance does not depend on the system size. There- 
fore, we can increase the radiation energy by adopting 
smaller samples with larger cavity frequency according 
to Eqs. 



In the present paper, only the fundamental frequency 
is kept for the plasma contribution to the total phase 
difference. This treatment is justified when the system is 
off resonance, where the amplitude of plasma oscillation 
is small and higher harmonics can be neglected safely. 

We have checked the validity of this treatment at the 
cavity resonance for the cavity mode (1,0) of the rect- 
angular mesa by comparing the estimates thus obtained 
with those by accurate simulations which include auto- 
matically all orders of higher harmonics Q]. We expect 
that this approximation also provides reasonable esti- 
mates for the (1,1) cavity modes of the rectangular mesa 
and for the cylinder geometry. 

In the resonating regime, the plasma part is much en- 
hanced, and higher harmonics appear [3j. The number of 
equations governing the kink state is 2m+2 when up to 
the m-th harmonics are covered. Here we provide equa- 
tions merely including the fundamental frequency with 
m = 1. While the formulas can be extended easily to 
higher harmonics, the resultant equations are much more 
complex, and take heavy numerics to solve. Discussions 
on this point will be reported elsewhere [14| . 

In the kink solution, the spatial part of the plasma 
term is the eigen function of the Laplace operator in the 
respective geometry with zero (or in presence of radiation 
very small) normal derivative at the edge, proportional 
to the magnetic field. In the cylinder geometry, it gives a 
cavity frequency voi/a with a the radius of cylinder and 
«oi = 3.3817 the first zero of J\. For a state presuming 
zero (or very small) electric field at the boundary, the 
typical frequency should be uoi/a, where uqi = 2.4048 is 
the first zero of Jo- Therefore, detecting the size depen- 
dence of the cavity frequency in the cylinder geometry 
can tell directly what state is realized inside the system. 
Since a rectangular mesa will give cavity frequency of 
ir/L for both of the two cases, the cylinder geometry is 
unique in determining the state realized in the junctions 
and thus the mechanism of the strong radiation observed 
recently in experiments. 

The present analysis reveals the way how the stack 
of Josephson junctions convert the dc energy to ac elec- 
tromagnetic radiation, by formulating out explicitly the 
oscillation amplitude and the frequency. The far-field ra- 
diation pattern is determined merely by the spatial part 
of the plasma term, and is the same as a thin capacitor 
subject to an ac driving voltage. Analyses on far-field 
radiation patterns for various modes are available in lit- 
erature [15| , and the results can be compared with those 
observed for the THz radiations from mesas of BSCCO 
single crystal. A detailed discussion will also be given in 
Ref. [3. 

To summarize, we have worked out explicitly the kink 
state in Josephson junctions of strong inductive coupling 
in 3D for both rectangular and cylindrical geometries. 
A set of equations are provided which permits one to 
understand the new state without heavy numerics. IV 



characteristics are revealed to be very nonlinear due to 
the cavity resonance where the plasma oscillation is much 
enhanced. The solution for a cylindrical mesa provides a 
higher resonating frequency than that of a square mesa 
with the same linear size by a factor of ~ 2.4. Investiga- 
tion on the size dependence of the resonating frequency 
in the cylinder geometry can give a direct evidence for 
the state realized in the stack of Josephson junctions, 
and thus can reveal the mechanism of the strong radia- 
tion discovered recently. Experiments using cylindrical 



mesas are highly anticipated. 
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